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Abstract 

We show that axinos, which are dominantly generated by the decay of the next-to-lightest super- 
symmetric particles produced from the leptonic Q-ball (L-ball) , become warm dark matter suitable 
for the solution of the missing satellite problem and the cusp problem. In addition, f2& — SIdm 
coincidence is naturally explained in this scenario. 
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I. INTRODUCTION 



Recent observations of cosmic microwave background anisotropies such as Wilkinson Mi- 
crowave Anisotropy Probe measured the abundance of components of the Universe very 
precisely. However, their origins are still one of the major mysteries of cosmology and parti- 
cle physics. The fact that the abundances of dark matter and baryon are of the same order 
may give us a great hint for their origins. 

In the minimal supersymmetric standard model (MSSM), flat directions consist of 
squarks/sleptons and produce a non-zero baryon or lepton number through the Affleck- 
Dine (AD) mechanism [lj]. Then, Q-balls, which are non-topological solitons 0], can be 
produced due to the instability and absorb almost all the produced baryon or lepton num- 
bers j3|. In the gravity mediated supersymmetry breaking model, the mass per charge of 
a Q-ball is larger than that of a nucleon so that they are unstable against the decay into 
light fermions. Then, they can directly decay into baryons and the lightest supersymmetric 
particle (LSP). In case that the charge of the produced Q-balls is large enough, Q-balls can 
survive even after the freeze-out of weakly interacting massive particles (WIMPs). Thus, 
the reason why energy densities of dark matter and baryon are almost the same magnitude 
can be explained in this scenario 0, lEf ■ 

However, it was pointed out that LSPs are often overproduced by the decay of Q-balls if 
the LSP is the lightest neutralino in the MSSM which gives the stringent conditions on 
the neutralino LSPs and AD fields. Only a few models free from this overproduction have 
been proposed @J8[. Instead, the supergravity models in which the LSPs is a stable gravitino 
are investigated [9]. In this scenario, Q-balls decay into the next-to-lightest supersymmetric 
particle (NLSP) directly instead of the LSP gravitino. The LSP gravitinos are produced 
by the decay of NLSP and becomes dominant over other gravitinos produced by thermal 
processed [lO, HI and by the decay of thermal produced NLSP H, 13[ in case that Q-balls 
can survive the evaporation. Then, it is found that such a gravitino dark matter scenario is 
still viable if the late decay of NLSP does not spoil the success of Big Bang Nucleosynthesis 
(BBN). Another interesting possibility is that the LSP is an axino, which is the fermionic 
superpartner of an axion. Axinos are also produced by the decay of NLSPs produced from 
the Q-ball [14 1 as well as by thermal processes and by the decay of thermally produced 

nlsp [H, mM, M, M ■ 



Such gravitinos and axinos often become an ideal candidate for cold dark matter. The 
models of cold dark matter (CDM) and dark energy combined with inflation-based scale- 
invariant primordial density fluctuations have succeeded at explaining many properties of 
the observed universe, especially the large scale structure of the universe. However, going 
into the smaller scales, some observations on galactic and subgalactic (< Mpc) seem to con- 
flict with predictions by high-resolution N-body simulations as well as analytic calculations 
based on the standard CDM model. The first discrepancy is called the missing satellite 
problem [2(j. The CDM-based models predict an order of magnitude higher number of 
halos than those actually observed within the Local Group. The other is called the cusp 
problem [2lj]. The CDM-based models also predict overly cuspy mass profile for the CDM 
halos compared to actual observations within the Local Group. In order to reconcile such 
discrepancies, several authors proposed modifications to the standard CDM-based model 



though the photoionization mechanism may overcome such difficulties [22J. One method is 



to reduce the small-scale power of primordial density fluctuations, which can be realized 



in a specific model of inflation [23|. Another is to change the properties of dark matter. 
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Spergel and Steinhardt introduced strong self-interaction among cold dark matter particles 
(collisional CDM), which enhances satellite destruction and suppress cusp formation [24]. 
The warm dark matter [25( , which can have relatively large velocity dispersion at the epoch 
of the matter-radiation equality can also reduce satellite production and cusp formation. 

In this paper, we consider axinos dominantly generated by the decay of NLSPs produced 
from the leptonic Q-ball (L-ball). Such axinos become warm dark matter suitable for the 
solution of the missing satellite problem and the cusp problem. In addition, Q b — VL DM 
coincidence is naturally explained through the Affleck-Dine mechanism and the subsequent 
L-ball formation in this scenario. In the next section, we discuss Q b — &>dm coincidence 
based on the Affleck-Dine mechanism and the subsequent L-ball formation. In section III, 
we show that axinos in our scenario become warm dark matter suitable for the solution of 
the missing satellite problem and the cusp problem. In the final section, we give concluding 
remarks. 



II. n b - n DM COINCIDENCE FROM AFFLECK-DINE LEPTOGENESIS 

We now discuss baryogenesis via Affleck-Dine leptogenesis and dark matter production 
from Q-ball decays within the framework of gravity mediated supersymmetry breaking. 



A. Lepton asymmetry 



The potential of the AD flat direction field is, in general, lifted by soft supersymmetric 
(SUSY) breaking terms and non-renormalizable terms @, 0. The full potential of the AD 
field is given by 

/J2~ 



v{4>) 



m 



Kin 



+ 



A 2 

(c 2 H + Am 3 / 2 ) A 



oiH' |0| 



nM 



n—3 



+ H.c. 



|2n-2 



+ A 2 



M 2n " 6 



(1) 



Here, is the soft SUSY breaking scalar mass for the AD field with radiative correction 
Kin |0| 2 . A flat direction dependent constant, K, takes values from —0.01 to —0.1 [28j]. A 
denotes a renormalization scale and —C\H 2 represents the negative mass squared induced by 
the SUSY breaking effect which comes from the energy density of the inflaton, with an order 
unity coefficient c% > 27|. A is the coupling of a nonrenormalizable term and M is some 
large scale acting as its cut-off. Terms proportional to A and c 2 are the A-terms coming 
from the low energy SUSY breaking and the inflaton-induced SUSY breaking, respectively, 
where m.3/2 denotes the gravitino mass. Here, we omitted possible terms which may appear 
by thermal effects 29|, |30j. These terms are negligible as long as we consider a sufficient 
low reheating temperature after inflation, as we will. Moreover, the model would face with 
"gravitino problem" 31], if the reheating temperature after inflation is so high that these 
thermal effect become effective, unless gravitino is LSP [§]. 

The charge number density for the AD field is given by n q = iq(<p*<f) — 0*0) where q is 
the baryonic (or leptonic) charge for the AD field. By use of the equation of motion of the 
AD field, the charge density can be rewritten as 



n q (t) 



a(ty 



dt'a(tr^^MA^ 



(2) 
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with a(t) being the scale factor. When the AD field starts to oscillate around the origin, the 
charge number density is induced by the relative phase between A-terms. By taking into 
account s = 4tt 2 g*T 3 /90, the charge to entropy ratio after reheating is estimated as 



T R n q 



4MpH 2 



q\A\Xm 3/2 T R \ 



n 



y os | 



M 2 p Hl s M n -z 



sin 5. (3) 



Here, Mp ~ 2.4 x 10 18 GeV is the reduced Planck mass, t os is the time of the start of 
the oscillation and sin<5 is the effective CP phase. In case that thermal corrections are 
ineffective, H os ~ m^, which yields 

^(^p 2 *. (4) 

From now on, as a concrete example, we consider a LLe direction of n = 6 as the AD field 
for our scenario. Since this is a pure leptonic direction, the lepton asymmetry generated by 
the Affleck-Dine mechanism can be estimated as 

^ ~ 1 x ip-M gW"** ( m ^ \ f ^y \ 3/2 / T R \ (M_\ 3/2 
s~ A 1 / 2 UoOGeW V J \100GeV J \M P J ' { ) 

B. Baryon asymmetry and LSP production from Q-balls 

The produced lepton asymmetry is not directly released to thermal bath. Instead, L-balls 
are formed due to the instability and almost all produced lepton numbers are absorbed into 
L-balls Q. 

First of all, we briefly summarize relevant properties of Q-balls in gravity mediated SUSY 
breaking models. The radius of a Q-ball, R, is estimated as R 2 ~ 2/(\K\mV) |5|. Numerical 
calculations provide a fitting formula for the Q-ball charge 



with 

e = — 



2q \A\^sm6 (7) 



where e c ~ KT 2 and (3 = 6 x 1(T 3 [3]. The Q-ball charge can be evaluated as 



(8) 



where we assumed e > e c because it looks to be more natural than the other which can be 
realized only for an accidental small sin 5 . Furthermore, if e < e c , additional "unnatural" 
parameters are required for our scenario, as we will show. 

A part of the charge of a Q-ball can evaporate by the interaction with particles in the 
thermal bath. The evaporation of charge of Q-ball is done by the evaporation with the rate 

Tevap = ~^ = -4vr J R Qj D cv n cq ~ -AnR 2 Q D cv fi Q T 2 , (9) 
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with D ev < 1 and by the diffusion with the rate 

r difr = -^ = -4nkR Q D difi n c * ~ -4irkR Q D diSf i Q T 2 , (10) 

where /xq is the chemical potential of Q-balls and the numerical constant k is very close to 
unity so that we will drop it hereafter. -D ^af ~ a/T is a diffusion constant [32j and a is a 
particle dependent coefficient given by jH, 34 

4 for squark 

6 for quark , . 

100 for left — handed (s)lepton 

380 for right — handed (s)lepton 

Here, we see that both the evaporation and the diffusion are efficient for low temperature, 
from Eqs. (jUJ) and ffTU]) with the relation between the cosmic time and the temperature: 

dT ivvfff for T<T " 

Moreover, by comparing r diff and r evap , 

~ ( m < . h W 1 \ (13) 



r \ T V 2 / V D 

1 evap \ ± I \ ev 



we can find that for low temperature 



1/2 



r<ay^~10 m ,(^)(|l) , (14) 

the diffusion is more crucial for estimation of the evaporated charge from Q-ball. Equation 
(flOl) is rewritten as 

§ = -^ZW* (|) . (15) 

Integrating Eq. ffT5l) from m^, because the evaporation from Q-ball is suppressed by the 
Boltzmann factor, we can estimate the total evaporated charge as 

AQ ~ 32kR Q J—a f j, Q ^^ 



3.2x2.4 ^/lTeVV/ T fl V / a \ /200 /0.01//i Q N . 

for ms>T R , 



AQ ~ x 10 * W (JL) ./M. /Ml f AfQ ) (17) 

/ 10 UJ ^300/ V 9* V 1*1 V™* 
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for < Tr. By taking Eq. ([8]) into account, we obtain 



Q 




(18) 



for < Tr 1 and find that about 10% of Q-ball charge would be evaporated. Here, one 
can see why the case of e < e c is irrelevant for us. If e < e c ~ 10~ 2 , e is replaced with e c 
in Eq. ffl8|) . Then these Q-balls cannot survive the evaporation unless the AD field mass is 
extremely small as m = O(10) GeV or (A 1/3 M P /M) 3 / 2 < 1 . 

The evaporated charges are released into the thermal bath so that a part of them is 



transformed into baryonic charges through the sphaleron effects [35(. Then, the resultant 
baryon asymmetry is given as 

n b 8 AQ n L 
~7 ~ 23 X ~Q X ~s 



-AQ x 10- 30 ( 15^XV' ( -^2—\ 
V ™><t> J V 100GeV / 



23 

\300l V 9, V 1*1 V«V \ (J^) form,, < T R 

Interestingly, the baryon asymmetry does not depend on the effective CP phase sin S unlike 
usual Affleck-Dine baryogenesis, because the CP phase dependences in both lepton asym- 
metry ni/s and the charge of Q-ball Q cancel each other. In addition, for < Tr, the 
baryon asymmetry basically depends on only one free parameter, the reheating temperature 
Tr, because other parameters are not free but known in a sense. 

When Q-balls decay, the supersymmetric particles are released from them. Since the 
Q-ball consists of scalar leptons, the number of the produced supersymmetric particles is 
given by 



Y ■ ■ N UL 
^NLSP — iv Q 



where Nq is the number of produced NLSP particles per one leptonic charge. 

Such produced NLSPs decay into axino LSP with a typical lifetime of 0(0.1 — 1) second. 
Thus, NLSPs produced by the Q-ball decay become a source of axino production. Of course, 
like gravitinos, axinos can also be produced by other processes such as thermal processes 
(TP), namely the scatterings and decays in the thermal bath, and non-thermal processes 
(NTP), say the late decay of NLSPs produced thermally. The relevant Boltzmann equations 
can be written as 

^NLSP + 3i?n N LSP = -(o^H^NLSP ~ n NLSp 2 ) + 7Q-baU ~ TnLSP^NLSP, (21) 

■ha + 3Hn a = (av(i + j -> a + ...))y-nin 3 - + (av(i -> a + ...))iTii + Tnlsp^nlsp, (22) 



1 For > Tr, the result is of the same magnitude but with different dependence on and Tr. 
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where 7Q_baii denotes the contribution to NLSP production by Q-balls decay, (crv)ij and 
(av)i are the scattering cross section and the decay rate for the thermal production of 
axinos, and Tnlsp is the decay rate of the NLSP. The total NLSP abundance, before its 
decay, is given by 

^nlsp = Nq— + ^nlsp, (23) 

where NqUl/s denotes the NLSP produced by L-ball decay and ^nlsp is the abundance of 
NLSP produced thermally and given by 



yTP _ H 
-^NLSP — „ 



T=m NLSP 



mNLSp/T/ . (24) 

\ I ann 



Here (crv) ann is the annihilation cross section and Tf ~ mNLSp/20 is the freeze-out temper- 
ature. The resultant total axino abundance is expressed as 

Yi = Y™ + Y™. (25) 

Here 

lf TP = W = Ar Q ^ + F N Y SP (26) 

is the nonthermally produced axino through the NLSP decay and denotes the axi- 
nos produced by thermal processes. For nonthermally produced axinos, while the NLSP 
abundance produced by L-ball decay is 

^ 2xl0 -,(^)(_^), (27) 

the typical value of l^xsp is given by 

NLSP V ^NLSP J V Mann / 

Thus, nonthermal production of axinos due to the thermal relic NLSPs decay, ^nlsp, can 
be negligible compared to that from Q-ball produced NLSPs, NqUl/s. On the other hand, 
axino production by thermal processes is dominated by scattering processes for the case 
that the reheating temperature is larger than the masses of neutralinos and gluinos. In this 
case, the abundance of such axinos is proportional to the reheating temperature Tr and the 
inverse square of Peccei-Quinn (PQ) scale f a and given by [l8| 



lTeV^ V fa/N J 



where N is the number of vacua and N = 1(6) for the KSVZ (DFSZ) model Thus, 



for Tr ~ 1 TeV, if f a /N > several xlO 11 GeV [38], F a TP is subdominant compared with 
1^ NTP ~ N Q n L /s. 

If this is the case, the energy density of axino is given by p a = m-a^NLSP due to the 
.R-parity conservation. Recalling 

PDA! _. ori .. n-10 ( ^DMh 2 



3.9 x 10- 10 ( ) GeV, (30) 
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the density parameter of axinos is expressed as 

Qah 2 ( ma \ ( n L / s \ ( N, 
0.11 



0.2GeW U x 10- 9 J 1 ' ^ 



Thus, axinos with the sub-GeV mass can be dark matter in our scenario. Now, one can see 
that the and Qdm is related through the lepton asymmetry. In fact, from Eqs. ( Tl9l) and 
( 13"Tj) . we obtain a relation between the abundances of dark matter and baryon asymmetry, 



Q b 2 fN Q \ ( AQ/Q \ ( 0.2 



tta 11 V 1 / V 1 x 10-1 / \m-a/m 



(32) 



where m p (~ 1 GeV) is the mass of proton. One may find the similar relation in the case 
of baryonic Q-ball (-B-ball) [13]. The difference between the case of .B-ball and L-ball is 
that the required mass of LSP from L-ball can be an order of magnitude smaller than that 
in .B-ball where the mass of LSP dark mater must be ~ 1 GeV, mainly because a part of 
lepton asymmetry produced by the Affleck-Dine mechanism, that is, only evaporated charges 
AQ / Q are converted to baryon asymmetry so that the number density of NLSPs produced 
by the L-ball decay become larger for a fixed baryon asymmetry. As shown in the next 
section, such a difference of axino masses is crucial for solving the missing satellite problem 
and the cusp problem. 

Equation f[3"Tj) with T R ~ 1 TeV to explain the observed baryon asymmetry yields a quite 
natural scale of the AD field mass 



™{l^)\£)W>- (33) 

As mentioned above, in our scenario, NLSP decays into axino at late time. Such late decay 
is potentially constrained by BBN. The lifetime of NLSP is given as 



for the case that the lightest neutralino x is NLSP in |18|]. Here, C a yY is the axion model 
dependent coupling coefficient between axion multiplet and U(l)y gauge field, Zu denotes 



the fraction of 6-ino component in the lightest neutralino. According to Ref. [181 
summerize the constraints as follows. First of all, for r x < 0.1 sec, there is no constraint. 
The corresponding mass of the NLSP neutralino is 

from Eq. ( 1341) . Thus, if m x > 320 GeV, this model is free from problems by the late decay 
of NLSP. This lower bound is a bit stringent than that in [18j], because we need to take the 
PQ scale somewhat larger as we mentioned. On the other hand, for 0.1 sec. < r x < 1 sec, 
the lower bound on axino mass exists and can be roughly expressed as 



> 

O.lGeV ~ 



(36) 
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by reading Fig. 4 in Ref . [18( . In this case, axino must be heavier than a few hundred MeV. 
For r x ~ 1 sec., the corresponding mass of the NLSP neutralino and the lower bound of 
axino mass is given by 

m x ~ 150GeV, 

m- a > 320MeV. (37) 



III. A SOLUTION TO THE MISSING SATELLITE PROBLEM AND THE CUSP 
PROBLEM 

An interesting consequence of such light axinos is their large velocity dispersion. There- 
fore, they can potentially solve the missing satellite problem and the cusp problem as stated 
in the introduction. For the scenario of dark matter particle produced by the late decay of 
long-lived particle, it is shown that the missing satellite problem and the cusp problem can 
be solved simultaneously if the lifetime of long-lived particle and the ratio of mass between 



dark matter particle and the mother particle satisfy the following relation [39 



6.3xl0 2 m a Y / ,/1.0xl0 3 m s \ 

sec < r x < sec, (38) 



m x J \ m x 



where we identified dark matter with axino (LSP) and the long-lived particle with the lightest 
neutralino (NLSP), respectively. 

Combining Eq. (1341) with Eq. (138!) . we have the following relation between the masses of 
the axino (LSP) and the lightest neutralino (NLSP), 



m ~ 



x , (0 ,3_0,3)(^)-'(^L_) (_^_)- 2 (_^_) 2 GeV . ,30, 



For rria — 1 GeV in the 5-ball case 14] . m x becomes a few GeV with its lifetime r x ^> 10 3 
second even if we take f a /N to be several times 10 11 GeV, which is excluded. On the other 
hand, for m„ = (9(0.1) GeV in the L-ball case of this paper, m x becomes (9(100) GeV with 
its lifetime r x < 1 second if we take f a /N to be several times 10 11 GeV. Thus, we find 
that axinos in this scenario can solve the missing satellite problem and the cusp problem 
simultaneously for natural mass scales with m„ = 0(0.1) GeV and m x = (9(100) GeV. 



IV. CONCLUDING REMARKS 

In this paper, we show that Affleck-Dine leptogenesis can explain baryon asymmetry 
and dark matter abundance simultaneously and that f2& — Qdm coincidence is explained for 
sub-GeV mass of the LSP axino. Though the basic idea is the same as Ref. (l4j, where 
S-balls are considered, the mass of LSP axino becomes an order of magnitude smaller in 
this scenario. On the other hand, the PQ scale is determined as f a /N = a few xlO 11 GeV, 
which will be tested in the future if PQ scale can be measured by e.g., the manner proposed 



m 



401 



The other attractive point is that axinos considered in this paper can potentially solve the 
missing satellite problem and the cusp problem simultaneously because they are relatively 
light and have large velocity dispersion. We have shown that axions in our scenario can solve 
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both problems for natural mass scales with m„ = 0(0.1) GeV and tbnlsp = 0(100) GeV. 
For simplicity, we concentrate on the case of the lightest neutralino NLSP with the mass 
to be (9(100) GeV. Such neutralinos are detectable in Large Hadron Collider. In addition, 
the corresponding lepton asymmetry is almost the maximal value under the assumption 
of (M 3 /\Mp) ~ 1. Hence, if this model is the simultaneous answer to both Q b — Qdm 
coincidence, the missing satellite problem, and the cusp problem, we would discover the 
NLSP neutralinos 2 with the mass of (9(100) GeV and their decays into photons (and 
negligible missing energy) with the lifetime of about (0.1 — 1) second. 
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